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A classical spin model derived ab initio for rare-earth-based permanent magnet compounds is presented. Our 
target compound, NdFe^N, is a material that goes beyond today’s champion magnet compound Nd 2 Fei 4 B in 
its intrinsic magnetic properties with a simpler crystal structure. Calculated temperature dependence of the 
magnetization and the anisotropy field agree with the latest experimental results in the leading order. Having 
put the realistic observables under our numerical control, we propose that engineering 5d-electron-mediated 
indirect exchange coupling between 4/-electrons in Nd and 3d-electrons from Fe would most critically help 
to enhance the material’s utility over the operation-temperature range. 

PACS numbers: 75.30.Gw, 75.50.Ww, 75.10.Hk, 75.10.Lp 


I. MOTIVATIONS 

In the past three decades, Nd 2 FeuBi has been the 
champion permanent magnet compound. A drawback 
of Nd 2 Fei 4 B has been its relatively low Curie temper¬ 
ature and that some practical applications require re¬ 
placements of Nd with heavy-rare-earth (HRE) elements 
such as Dy and Tb to enhance the high-temperature co- 
ercivity, which is roughly proportional to the anisotropy 
fields. Since the HRE elements are less abundant, ways 
to achieve the equivalent magnetic properties to those 
of HRE-doped permanent magnets using only light rare- 
earth (LRE) elements have recently been sought after—. 
Also the relevance of understanding and controlling the 
finite-temperature magnetism of 4/-3d intermetallics is 
appreciated ever more. 

A possible solution was recently suggested 1 by the 
material NdFe^N stabilized in an almost bulk stated, 
where a newly fabricated film sample that consists of 
more than hundred unit-cell layers shows superior intrin¬ 
sic magnetization and anisotropy field to Nd 2 Fei 4 B. The 
materials family, RFei 2 - x T x (N) (R=rare earth), had ac¬ 
tually been known for a long time^ where the nitrogena- 
tion pulls up the Curie temperature by 100-200 [K] and 
the magnetic anisotropy is enhanced as well, but the 
achieved magnetic properties were not on a par with the 
champion magnet compound Nd 2 Fei 4 B partially due to 
the necessity for the presence of the third element T=Ti 
etc. to stabilize the particular crystal structure. Recent 
breakthrough^ made it possible to have NdFei 2 without 
the third element in a sample fabricated as a thick film 
and nitriding achieved the intrinsic magnetic properties 
that goes beyond Nd 2 Fei 4 B at high temperatures^. 

Thus we are motivated to theoretically address the 
finite-temperature magnetism of NdFe^N. This would 
provide the prospect for its intrinsic magnetic proper¬ 
ties, which serves to solve the high-temperature coer- 


civity problem in LRE-based compounds. On the the¬ 
ory side, finite-temperature magnetism of rare-earth- 
based permanent-magnet materials poses a fundamen¬ 
tally challenging many-body problem: ab initio predic¬ 
tions mostly focus on the ground-state properties and 
the finite-temperature magnetism was discussed at best 
on the basis of a mean-field theory of a simplified model 
on the basis of a molecular field acting on an isolated 
rare-earth magnetic moment—. For comparison with 
experimentally observed magnetic anisotropy, contribu¬ 
tions from 3d electrons are sometimes added in an ad 
hoc manner. In principle, the theory of magnetism in 
4/-3d intermetallics takes a description of the correlated 
electrons in 4/ and 3d-orbitals, which may be done with 
a multi-orbital periodic Anderson model (PAM) with the 
conduction electrons composed of 5d-band out of the 
rare-earth elements, harboring two species of impurities, 
3d and 4/, each with the different levels of on-site elec¬ 
tronic correlation. 

In order to meet the urgent practical needs and also 
to provide a guideline data for future realistic simula¬ 
tions of PAM with huge number of orbitals, we exploit 
the essence of a simplified models to describe only the 
low-energy effective physics of 4/-3d intermetallics with 
the model parameters determined as realistically as pos¬ 
sible through ab initio calculations: we define a multi¬ 
sublattice spin model with one group of the sublattices 
describing the 4/-originated localized magnetic moments 
and the other describing the 3d-magnetization; the spins 
reside on a realistic lattice that mimicks the crystal struc¬ 
ture of the given target material NdFe^N with the ab 
initio input parameters. Note that the effective parame¬ 
ters, such as the strength of exchange couplings and crys¬ 
talline electric fields, are determined as a consequence of 
the electronic states. They are target-material depen¬ 
dent, thus first-principles evaluation of the parameters 
is crucial for the quantitative modeling. Then we solve 
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the finite-temperature many-body problem with numer¬ 
ically exact Monte Carlo method to get the temperature 
dependence of magnetic observables and quantitatively 
compare with the latest experimental data. Our realis¬ 
tic lattice model is more realistic than were discussed in 
the previous works: The sublattice-specific character of 
each of the Nd and Fe atoms in the unit cell is taken 
into account on the basis of first principles, which is in 
contrast to a uniform molecular field imposed by Fe act¬ 
ing on rare-earth magnetic moments. Establishing the 
computational control of the intrinsic properties of mag¬ 
netism of NdFei 2 N, we discuss within our model how to 
manipulate it to enhance its practical utility most effec¬ 
tively. 

II. THE REALISTIC-LATTICE SPIN MODEL 

The spin model Hamiltonian defined on the lattice of 
the given crystal structure reads as follows. 


H = Ht + "Hr + 'Hrt- 

(1) 
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We have denoted the lattice points on which 3d (4/) mag¬ 
netic moments reside by T (R), respectively. Here S,; is 
a magnetic moment defined on site i, Jij is the exchange 
coupling between localized magnetic moments S 7 ; and S j, 
and Di encodes the single-ion magnetic anisotropy energy 
(MAE). Note that the summation runs over each 

bond connecting the sites i and j only once. The 4/-part 
is described explicitly with the total magnetic moment 
J = L -f S with L being the orbital moment. The spin 
moment can be extracted via S = ( gj — 1) J with gj being 
Lande’s g-factor. Describing a 4/-mom ent of Nd 3+ with 
a classical spin of length gj\f J(J + 1) with J = 9/2 is 
semi-quantitatively justified within the scope of setting 
the target temperature range in 200 [K] ~ T ~ 500 [K], 
with the J-multiplets separated in the scale of 1000 [K]. 
The 4/-4/ exchange coupling terms, J RR , in the scale of 
0(1) [K] have been dropped because our target energy 
scale to be realistically described is motivated by the typ¬ 
ical operating temperature range 200 [K] ~ T ~ 500 [K] 
of permanent magnets. Magnetic energy scales coming 
from J RR ’s are at most ~ 0.1T and should be washed 
out. 

The 4/-3d indirect exchange coupling as denoted by 
J R J in Eq. ©) comes from (5d) m -(3d) n ( m ~ 1 and 
n > 5) exchange which is antiferromagnetic and RE on¬ 
site 4/-5d direct exchange which is ferromagnetic. Over¬ 
all J R J is an antiferromagnetic coupling between the spin 
component of 4/ and 3d, which means 4/ total moment 
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FIG. 1. Crystal structure of NdFei 2 N. Along the Nd(2a)- 
N(26)-Nd(2a) direction runs the c-axis. Perpendicular to it, 
the other two equivalent directions, a-axis and b-axis, span 
the Nd(2a)-Fe(8i)-Fe(8i)-Nd(2a) lines. 

and 3d magnetization are ferromagnetically coupled for 
Nd 3+ with gj = 8/11. The overall scale factor o:rt, 
tentatively set to be one, has been introduced to phe¬ 
nomenologically describe the indirect nature of the 4/-3d 
exchange coupling. 

A. Derivation of the leading-order parameters 

As shown in Fig. [L] the material NdFe^N has a 
ThMni 2 crystal structure (space group M/mmm ) with 
Nd occupying the body-centered site of a tetragonal unit 
which incorporate 2 formula units. Around each Nd site, 
there are 4 Fe(8i) atoms which surround along the a- 
axis and b-axis. Four Fe(8j) atoms make a square right 
in the middle of Nd sites along the c-axis and nitrogen 
goes in the center of this square. Fe(8/) sites make a 
tetragonal box with the center occupied by Nd. Tak¬ 
ing the ThMni 2 crystal structure as the starting point, 
the crystal structure is optimized from first principles. 
Then the crystal field parameters on Nd 3+ are calcu¬ 
lated on the basis of open-core description for the 4/- 
shelli^. Up to here the calculations are done with the 
ab initio electronic-structure computational code pack¬ 
age, QMASi£. The parameters for the exchange-coupling 
are calculated with the ab initio electronic-structure cal¬ 
culation code, Machikaneyama (AkaiKKR)i£, using the 
Korringa-Kohn-Rostoker Green’s function method, fol¬ 
lowing the prescription proposed by Liechtenstein et 
al.—. Re-writing the term in Eq. © as J^ T S, • S j = 
(SiJ^ T Sj)e i ■ ej with the vector e 7 denoting the direc¬ 
tion of the magnetic moment on site i, calculated ex¬ 
change couplings Sj) are summarized in Fig. [2] (a) 

as a function of inter-site distance. As is the case for the 
champion magnet compound Nd 2 Fei 4 Bi 2 , we find that 
the dominant Fe-Fe exchange couplings come from Fe 
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FIG. 2. (a) Calculated Fe-Fe exchange couplings as a function 
of Fe-Fe separation for NdFei 2 and NdFei 2 N. (b) The same 
data set is plotted on a plane spanned by Jij (NdFei 2 ) and 
Jij (NdFei 2 N), with Jij shared among NdFei 2 and NdFei 2 N 
as the spatially equivalent bonds. 


pairs whose interatomic distances almost coincide with 
that in ct-Fe. In our model we incorporate those domi¬ 
nant Fe-Fe exchange couplings, that can be classified into 
8 classes on the realistic lattice for NdFei 2 (N) to make 
a leading-order description. Comparing the calculated 
exchange couplings on the same bonds located on the 
equivalent positions in the unit cell between NdFei 2 and 
NdFei 2 N as shown in Fig. [2] (b), nitrogenation-enhanced 
exchange couplings are identified. Interestingly, we ob¬ 
serve those bonds are located along a kagome-lattice net¬ 
work spanned by Fe(8/) and Fe(8j), reminiscent of one- 
generation-back permanent-magnet materials family rep¬ 
resented by SmCos with the strong magnetic anisotropy. 
The other leading-order term in the overall Hamiltonian 
written as Eq. © is 4/ single-ion MAE D R (J^) 2 , which 
also shows up in the energy scale of 0(10) [meV]. Re¬ 
writing the term in Eq. m as 

Dl {J z m ) 2 = Kl (e m • n ) 2 , (5) 

with n denoting the direction of the easy axis and e the 
direction of the total magnetic moment J m on the rare- 
earth site m, a QMAS calculation gives AT Nd = 5.8 or 
11.0 [meV] (easy-axis) for NdFe^N and A' Nd = —3.0 or 
—2.3 [meV] (easy-plane) for NdFei 2 depending on the 


TABLE I. Calculated 4/-3d exchange couplings in meV. 
The notation follows Eq. ©■ The site index m denotes the 
Nd(2a) sublattice. 


i 

NdFei 2 NdFei 2 N 

Fe(8i) 

1.58 

1.90 

Fe(8j) 

1.54 

1.00 

Fe(8/) 

1.93 

2.36 


calculation setups^. As a first step, we choose I\ Nd = 
8.1 [meV] for NdFe^N and = —2.8 [meV] for NdFei 2 
within the range. 


B. Sub-leading parameters 

We define J RT and D T in the overall Hamiltonian writ¬ 
ten as Eq. ©■ Calculated antiferromagnetic exchange 
coupling between 5d-bands of Nd and 3d-bands of Fe 
are taken as they are to be the 5d-mediated 4/-3d cou¬ 
pling, assuming that intra-atomic direct exchange cou¬ 
pling is big enough to let the localized 4/-moment and 5d- 
polarization work as a unified body. The results are sum¬ 
marized in Table |U The Fe-originated MAE in Eq. ([2]), 
which we will denote analogously to Nd-originated one 
in Eq. m as follows, 

( 6 ) 

comes in the order 0(0.1) [meV] as referring to the past 
experimental measurements^ done for YFenTi. In the 
present modeling we just incorporate K Fe = 0.1 [meV] 
uniformly on all Fe sites as a phenomenological setting 
which should be sufficient to describe the magnetic prop¬ 
erties around the room temperatures and slightly higher. 

Higher-order terms in addition to AT Nd in Eq. m can 
be suspected in principle for which we have also calcu¬ 
lated and saw that they are an order of magnitude smaller 
than the leading-order ones. Within the present scope to 
pick up the leading-order behavior of finite-temperature 
observables, we have dropped the higher-order contribu¬ 
tions to the single-ion magnetic anisotropy of Nd. 


C. Methods 

Having defined the realistic model basically from first 
principles for NdFei 2 (N), the temperature dependence of 
the magnetic properties are calculated using the classical 
Monte Carlo method with the Metropolis local updates. 
We do one of the most plain local updates, that is, picking 
up lattice sites stochastically and proposing uniformly on 
the spherical spin space up to the stochastic decision re¬ 
ferring to the energy difference as calculated by the spin 
model Hamiltonian defined as Eq. ©■ Effective one lat¬ 
tice sweep, i.e. with the stochastic choice of the lattice 
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site to be updated done N S { te times, is counted as one 
Monte Carlo step. Our input parameters are simply set 
to the T = 0 ab initio values. For typical runs in the 
present work, 10 5 Monte Carlo steps (MCS) with each 
step consisting of a lattice sweep were sufficient to reach 
the thermal equilibrium state and subsequent 10 6 MCS 
are used for the numerical measurement of observables. 
The tetragonal unit of the body-centered-network of Nd 
is counted as a system-size unit in our calculation, thus 
a calculation with L x L x L system actually contains 
26L 3 = iV s i te spins with the tetragonal unit having two 
formula units. Majority of the data are taken with L = 4 
which is sufficient as seen below with Fig. [4] (a) while 
the data points on the high-temperature side took sys¬ 
tem sizes up to L = 10, with the anisotropy field getting 
smaller and the calculated system getting closer to the 
Curie temperature where the finite system size can be¬ 
come an issue due to the diverging correlation length. 


III. RESULTS 

With the prescription for the construction of the 
realistic-lattice spin model, calculated temperature de¬ 
pendence of magnetization and anisotropy field are 
shown in comparison with the experimental data£. Model 
parameter dependence of the overall calculated temper¬ 
ature dependence is inspected to draw our main conclu¬ 
sion: 4/-3d exchange couplings in Eq. JI|) dominates 
the observables in the operation temperature range i.e. 
around the room temperature or higher. 


A. Magnetization 

Calculated temperature dependence of magnetization 
M(T ) = ^M x 2 + My + M? for the bulk, the Nd- 

sublattice, and the Fe-sublattice of NdFe^N is shown in 
Fig. El (a). Calculated Curie temperature falls in the same 
range as the experimentally claimed one. Considering the 
possible presence of a-Fe-originated noise in the experi¬ 
mental data and neglect of all next-nearest-neighbor ex¬ 
change interactions in our model which should have low¬ 
ered the computational magnetic-ordering energy scale, 
the agreement is satisfactory. Also we see that the 
system-size dependence is negligible in our focus tem¬ 
perature range, which up to T = 500 [K] at most. 
Comparing the calculated temperature dependence of 
magnetization between NdFei 2 and NdFe^N as shown 
in Fig. [3] (b), we see that the experimentally observed^ 
nitrogenation-triggered enhancement in the Curie tem¬ 
perature by ~ 200 [K] is well reproduced. 

We note that Nd in NdFei 2 has the easy-plane mag¬ 
netic anisotropy^i^ which would compete against the 
easy-axis anisotropy from at least a part of the Fe sublat¬ 
tices. Here we just track the origin of the difference in the 
Curie temperature between NdFei 2 and NdFe^N to the 


(a) 



T [K] 


(b) 



T [K] 

FIG. 3. (a) Calculated temperature dependence of magnetiza¬ 
tion for NdFei 2 N compared to the latest experimental data£. 
(b) Comparison of the calculated temperature dependence of 
magnetization between NdFei 2 and NdFei 2 N within L — 4. 


nitrogenation-enhanced exchange couplings as demon¬ 
strated in Fig. [2] Further studies on the outcome of the 
competing anisotropies in NdFei 2 is separated for future 
work. 


B. Anisotropy field 

Calculated magnetization curves with the externally 
applied magnetic field parallel and perpendicular to the 
easy axis in NdFe^N (which is c-axis) at T = 348 [K] is 
shown in Fig. [2(a). The anisotropy field H a is identified 
as a crossing point of the two curves with a linear fit. 
Our numerical measurement just follows the experimen¬ 
tal way to determine H a . Thus determined anisotropy 
field H a (L,T) with the L-dependence being saturated 
out within the statistical error bars is plotted as a func¬ 
tion of temperature in Fig. [4] (b) in comparison to the 
recent experimental data^. In the operation tempera¬ 
ture range, a leading-order numerical control seems to 
be achieved. 
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FIG. 4. (a) Numerical determination of the anisotropy field 
via the calculated magnetization curves for NdFei 2 N at T = 
348 [K] and (b) thus calculated temperature dependence of 
the anisotropy field for NdFei 2 N. 


C. Roles of each parameter 

We inspect which parameter dominates which part 
of the temperature dependence of the observables com¬ 
paring with the latest experimental data. One of the 
most important practical utilities of permanent magnets 
lies in the coercivity that is roughly proportional to the 
anisotropy field of the parent material or the main phase, 
and the main issue here is to figure out how to sustain the 
anisotropy field at high temperatures. Calculated MAE 
for Nd, A' Nd , and the input parameters for K Fe obviously 
influence the temperature dependence of the anisotropy 
field. Flipping the sign of K Fc for the Fe sublattice to 
have an easy-plane anisotropy, the overall reduction of 
the anisotropy field by 2 x 12 x |AT Fe |/A' Nd ~ 30% at 
the lowest temperature range is observed as shown in 
Fig. [5] (a). An interesting observation is made by en¬ 
hancing AT Nd by 1/3 and flipping the sign of AT Fe at the 
same time, where the overall anisotropy field at the low- 
temperature side almost collapses onto the original tem¬ 
perature dependence as shown in Fig.[5](b). The contri¬ 
bution of K Fe is found to be quantitatively important in 
the bulk MAE at finite temperatures where an increase 
in the anisotropy field only by 10% can lead to a practi¬ 
cal breakthrough in designing LRE-based magnets^. Ab 


( a ) 



(b) 



T [K] 

FIG. 5. (A Nd , A' Fe )-trends of the calculated temperature de¬ 
pendence of the anisotropy field for NdFei 2 N. (a) Trends of 
H a (T) with respect to K Fe = ±0.1 [meV]. Data with K Fe = 0 
is also included as a reference, (b) A Nd is enhanced by a frac¬ 
tion of 1/3 in observing the trends of H a (T) with respect to 
K Fe = ±0.1 [meV]. 

initio determination of the exact nature of Fe-sublattice- 
dependent MAE is under way. 

Numerically manipulating the 5d-mediated 4/-3d indi¬ 
rect exchange coupling, J RT , we observe how the calcu¬ 
lated temperature dependence of the magnetization and 
the anisotropy field in NdFe^N is affected as shown in 
Figs. [G] (a) and (b), respectively. The anisotropy field is 
more significantly affected by J RT than the magnetiza¬ 
tion is especially around T ~ 300 [K]. This is a numerical 
demonstration that a small enhancement in J RT , which 
may be realized by the conduction-band engineering in 
4/-3d intermetallics, can already lead to a considerable 
improvement. 


IV. DISCUSSIONS 

Having identified the quantitative relevance of J R J 
written as Eq. © for the magnetic properties of NdFe^N 
in the operation temperature range, we further pin-down 
the role of in determining the realistic magnetic 
anisotropy energy scale at finite temperatures, through 
a close comparison between calculated results and the 
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(a) 



T [K] 


(b) 



FIG. 6. J RT -dependence of the calculated temperature de¬ 
pendence of (a) magnetization and (b) anisotropy field for 
NdFei 2 N. Given values from the ab initio electronic-structure 
calculations are tabulated in Table [T] 


experimental data. 

In Fig. E it is seen that the experimentally observed 
slope and the downward convexity of the temperature de¬ 
pendence of the anisotropy field for NdFe^N is closest to 
the calculated H a (T) with Ort = 0.25. Also the experi¬ 
mental anisotropy field near the liquid-nitrogen tempera¬ 
ture seems to come close to 14 [T]2i which suggests possi¬ 
ble systematic underestimate in our calculated magnetic 
anisotropy which falls below 10 [T] in the limit T —> 0. 
The ab initio estimation for the crystal-field coefficients 
for the estimation of the uni-axial magnetic anisotropy 
energy indeed involves a certain uncertainty^. 

An improved data collapse between theory and exper¬ 
iment in the lowest temperature range of H a (T ) can be 
observed by manually setting AT Nd = 2A// d . As shown in 
Fig.jTJ the slope of the temperature dependence up to the 
room-temperature range is well reproduced. The upper 
shift of the calculated H a (T) on the high-temperature 
side of the operation temperature range can be adjusted 
by a manual scaling of J TT with an overall factor of 4/3 
to match the Curie temperature. This can be considered 
as an effective renormalization of the 3d-3d exchange cou¬ 
plings imposed by the discarded longer-range part in the 
exchange couplings. 

Thus an inspection of the experimental temperature 



T [K] 

FIG. 7. For NdFei 2 N, a data set which describes the experi¬ 
mental data most quantitatively comes with «rt = 0.25. 

dependence of magnetization and anisotropy field leads 
to a set of model parameters for a quantitative descrip¬ 
tion of the experimental finite-temperature data. Based 
on the obvious relevance of J TT for the Curie tempera¬ 
ture and AT R for H a (T = 0), it is seen that J RT deter¬ 
mines the slope of H a (T) near T ~ 0. 

In such parameter set, the 4/-3d exchange couplings 
come close to a quarter of the tabulated numbers which 
were 1 ~ 2 [meV] as seen in Table HI In the language of 
the two-sublattice models where an isolated 4/-electron 
magnetic moment is put into the sea of 3d-electron mag¬ 
netization via a molecular field H m , the molecular field 
can be written in terms of our lattice model language as 
follows. 

= aRT2J RT S x ,z Nd 

where ZNd is the coordination number around the Nd 
magnetic moment which is in the present case ZNd = 20 
within the nearest neighbor. Plugging in the realistic 
numbers anr ~ 0.25, J RT ~ 2 [meV], S ~ 2 [/tb] the 
magnitude of the exchange field is found to be H m ~ 
40 [meV] = 700 [T] which gives the same order as was 
found in the experimental analyses H m ~ 450—600 [T] by 
neutrons^. <trt = 1 corresponds to H m beyond 1000 [T] 
thus is out of the realistic scale. The smallness of the 
realistic number in the factor curt presumably reflects 
the indirect nature of the 4/-3d exchange. 

V. CONCLUSIONS AND OUTLOOK 

Ab initio modeling for the 4/ and 3d magnetism cou¬ 
pled by 5d-electrons on the basis of the realistic spin- 
lattice model for the rare-earth permanent magnet ma¬ 
terials NdFei 2 N and NdFei 2 quantitatively captures the 
realistic energy scales in the leading order in the oper¬ 
ating temperature range, 200 [K] ~ T ~ 500 [K]. Ex¬ 
perimentally observed magnetic ordering and magnetic 
anisotropy energy scales are put under numerical control 
and we have shown that enhancing the 4/-3d indirect 
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exchange coupling would work most effectively to realize 
the magnetic properties of more practical use in the op¬ 
eration temperature range. Establishing a quantitative 
description starting with the simplified model would pave 
the way to more realistic simulations to explicitly incor¬ 
porate the strongly-correlated nature of 4/-electrons em¬ 
bedded in the conduction-electron sea of 5d-electrons to 
comprehensively describe the finite-temperature physics 
starting from the range T ~ 1 [K] all the way to the 
Curie temperature close to 1000 [K]. 
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